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The dynamics of the quantum entanglement is a fundamental characteristic for various quantum 
systems. Since the computable entanglement measure for higher dimensional quantum states itself 
is absent, the dynamics of the entanglement expressed in an operational method will be of interest. 
We study the dynamics of r, an analytical lower bound of squared concurrence, of a bipartite d® d 
quantum state when one party goes through an arbitrary noisy channel. For a pure input state, the 
range of r is obtained explicitly. For a mixed input state, an upper bound of r is found. Interestingly, 
the tangle r', as an upper bound of squared concurrence, also has a similar dynamical property. Our 
results are similar to that of Konrad et al. and can help the estimation of high-dimension bipartite 
entanglement in experiments. 
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I. INTRODUCTION 

Quantum entanglement, which is considered to be the 
most non-classical phenomenon in the quantum world, 
lies in the central position of quantum information the- 
ory (QIT). It has been identified as a key resource in 
many aspects of QIT, such as quantum teleportation, 
quantum key distribution and quantum computation [l[ . 
But while implementing quantum information process- 
ing in real physical systems, it's inevitable that the en- 
tanglement decays due to the interactions of our system 
with the environment, making it significantly important 
to study the dynamical property of entanglement in re- 
alistic situations. 

The dynamical property, namely the time evolution of 
entanglement of a state is usually deduced from the time 
evolution of the state itself [2] 3]. However, recently, in 
Ref.fl], without solving the master equation of a quan- 
tum state but by utilizing the Jamiolkowski isomorphism, 
Konrad et al. presented a factorization law for a two 
qubit system, which describes the evolution of entangle- 
ment in a simple and general way. Then, Li et al. gener- 
alized this result to that of a bipartite quantum system 
of arbitrary dimension ||. In the study above, concur- 
rence which is a well accepted entanglement measure, was 
used to quantify the entanglement. As is well known, for 
higher dimensional bipartite quantum state, there is no 
analytic method in general to find concurrence. Thus it 
will be very interesting if we can study the dynamical of 
the entanglement which is quantified in an operational 
way. Unfortunately, there is no such an operational mea- 
sure of entanglement for an arbitrary bipartite quantum 
state. However, there exist a lower bound of squared 
concurrence which is analytic [|| and we represent it as 
t. In this paper, we will investigate the dynamical of the 
lower bound of the squared concurrence t. As a special 



case for two-qubit state, our result reduces to the result 
by Konrad et al. Moreover, the tangle r' as defined in 
Rcf. [7] is an upper bound of squared concurrence. Inter- 
estingly, it has a similar dynamical property with r. To 
clarify our results, we use the depolarizing and the phase 
damping channels as the examples. 



II. CONCURRENCE AND ITS UPPER AND 
LOWER BOUNDS 

As the beginning we recall the definition of con- 
currence, t and tangle t' . For a pure bipartite 
state pab —\ ip)(ip\ in a finite d\ ® g?2 dimensional 
Hilbert space " Ha ® H-b, t he concurrence is defined 
as C(|V)) = v 72 ^ -Trp 2 A ), with p A = Tr B p A B the 
reduced density matrix. For a mixed bipartite state 
P = 52iPi\ipi){tl>i\,Pi > 0, Y^iPi = 1) the concurrence 
is defined as the convex roof of all possible decompo- 
sitions of p into the pure states namely C(p) = 

™ a -{p,M.)}Y Ji P l C{\'il ) *i))- 

Although the concurrence of a general bipartite mixed 
state defined above is difficult to solve due to a high- 
dimensional optimization, a computable lower bound of 
squared concurrence can be found in Ref.[6]: 



C 2 (P) > 



di-l d.2-1 

Y Y C pr,p'r'(p) = t (p), 
r>p=0 r'>p'=Q 



(1) 
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where r is a lower bound of squared concurrence and 

Cpr,p'r'(p) max{0, ^p r p' r r ^pr.p'r' ^pr.p'r' ^pr,p'r' }' 

(2) 

with X pr , r i being the squared roots of the four nonzero 
eigenvalues, in decreasing order, of pp pr ,p' r ', where 
Ppr, P 'r' = (L pr ® L p > r >)p*(L pr <gi L p i r >) and L pr = \p)(r \ - 
\r)(p\(p,r = 0,l,...,di - l;p < r), L p , r , = \p')(r'\ - 
\r'}(p'\(p' ,r' — 0,1, ...,d 2 — l;p' < r') are the genera- 
tors of the group SO(di) and SO(d2) respectively. It's 
clear that r can always be calculated analytically. 
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According to Ref. [6], every C prtP i r i(p) can be seen as 
a two qubit concurrence of a 4 x 4 matrix p, which is a 
submatrix of p, 

(Ppp' ,pp' Ppp 1 ^pr' Ppp' 7 rp f Ppp' ,rr' \ 

Ppr' ,pp' Ppr' ,pr' Ppr' ^rp' Ppr' ,rr' | 

Prp' ^pp' Prp' ,pr' Prp' ^rp' Prp' ,rr' I 

Prr' ^pp' prr' ^pr' Prr' ,rp' Prr' ,rr' / 

with subindiccs p and r associated with TIa and p' and 
r' with Hb- So r(p) in fact is the sum of some two qubit 
entanglement in a high dimensional state, according to 
which we can rewrite Eq.(l) in another form: 



v 



r(p)^^C 2 (p 4 ), 



(4) 



where C is just the two qubit concurrence and pi is a 
submatrix of p of the form (3), the number of which is 

-T-) _ did 2 (di-l)(d 2 -l) 
u — 4 

One can prove that t is a convex function of the den- 
sity operator. According to the definition, Ei-Pi/ 9 *) = 

EfcLi C 2 ^ J2i PiPi^j > where p\ is a submatrix of p; t . Using 

the convexity of C, T^iPiPi) < ELi (l>2iPi C (Pi)) ■ 
Recall that f(x) — x 2 is a convex function, namely 

V 



k=l 



-(E^)^EE^ 2 ^) 

i 

^p l r{p i ), (5) 



which is just what we want to prove. 

The tangle for a general mixed state is defined as 

/(p)= min J>C 2 (hM), (6) 
{p., I'M} ^ 

which is also a convex function of the density operator. 
One can easily see that r'{p) > C 2 {p) from the convexity 
of concurrence so it's an upper bound of squared concur- 
rence. 



III. THE DYNAMICS OF CONCURRENCE 

Here we briefly review the dynamics of concurrence 
demonstrated in Ref. [4] and Ref. [5]. First consider a 
2 ® 2 two qubit pure state, after only one qubit goes 
through an arbitrary channel £, the concurrence be- 
tween them decays just by a universal factor only de- 
termined by £'s action on the maximally entangled state 
l^ + ) = 75(100) + |11)): 

c[(i ® emm = c[(i ® £)w + ){<t> + \]cm). (?) 



For a general d\ ® di pure state, a similar relation is 
satisfied, with a sacrifice that the equality is replaced by 
an inequality: 



C[(1®£MM] < f C[(l®5)|^+)<0+|]C(^». (8) 



Both results above can be generalized to the case where 
the input state is a mixed state. For a 2 ® 2 mixed state 
we have 

C[(l ® 5)po] < C[(l ® £)\4> + )^ + \]C(p») (9) 
and for a d\ ® d 2 mixed state we have 



C[(l ® < yC[(l ® £)|0 + )(^+|]C(po). (10) 



IV. THE DYNAMICS OF r AND r' 



Generally speaking, to solve the concurrence of a high- 
dimensional mixed state, just like (1 ® £)\4> + )(<fi + \, one 
must make an optimal decomposition of the state, which 
is a notoriously difficult task, making the right hand side 
of Eq.(8) and Eq.(10) nearly impossible to be calculated 
analytically except in some special cases. This motivates 
us to investigate the time evolution of r, which can be 
calculated analytically. 

Let us consider a d<g>d bipartite quantum system whose 
Hilbcrt space is 7i, then any pure state \ip) £ H can be 
expressed by Schmidt decomposition as follows: 



d-l 



d-1 



iv>) = E i^Z^i = 1 



(ii) 



and the maximally entangled state in Tl can be written 

Because Jamiolkowski isomorphism can be extended 
to bipartite systems of arbitrary finite dimension, the 
dual picture used in Ref. [4] will be valid in higher di- 
mensions. Consider a quantum channel £, according to 
Jamiolkowski isomorphism, when only one qubit of the 
state (11) goes through £, we have p' = = 

where p £ = ag£H ^£l with p, p> and p" 
the normalization coefficients and one can verify that 
p' = d 2 pp" . The action of channel £^ can be expressed 
in a simple form that ® l)pe = (A4 ® \)pe{M^ ® 1), 

where M = ^ Eto V<*l*}<*l- 

Because = M and ML pr M = ^ UJ ? UJr L pr , 

&ct(p' p' prp , r , - Al) = dot {^rp £ P£ pr , p , T , - Al), from 
which we can have that C prp , r ,((l ® £)\ip){i})\) = 
d 2 u p uj r C 2 rip , r ,((l ® £)\<f> + ){<P + \). Noting C 2 prip , r ,m = 
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AuipOj r 5 pp i5 rr i , an important relation is derived: 

r">p"=0 

xC P ». r »,pv((l®f)|0 + )<0 + |)) 2 .(12) 

In Introduction we have explained C pr y r > as a two 
qubit concurrence, so Eq.(12) means that the evolution 
of a certain two qubit entanglement in a high dimen- 
sional state also obeys a law which is similar to Eq.(7) 
but more complicated because we must consider all re- 
lated two qubit entanglement, as demonstrated in the 
sum in the RHS of Eq.(12). It's easy to see that when 
d = 2, Eq.(12) is equivalent to Eq.(7). 

In what follows we want to find the range of r((l £g> 
£)\ip)(ip\). According to the definition of r, we have 

t((i®£M(vd = ^c 2 (iv>)) e E ^ d T r 

r>p=0r<>p<=0 l^i<j=0 ^i^j 

xC 2 pr , p , r ,((l^£)\cj ) +)(^\). 

(13) 

Considering u> p uj r < J2i<j=o UJiUJ j' we immediately get 
the upper bound of r((l ® £ )\ip)(i/)\): 

r((l®£)|^)(^|)<^r((l®£)|0+>(0+|)C 2 (|^)). (14) 

On the other hand, one can show that Yli<j=o ^i^j = 

K 1 ~ J2i=o u i) ^ ^T- Let ? 7 = min {p,r} u p u> r for any 
pair p < r satisfying ui p uj r ^ 0, we find a lower bound of 

t((1<8)£MM): 

r((l®£)|^)(V|)>|^^r((l®£)|0 + )(0 + |)C 2 (|^). 

(15) 

Eq.(14) and Eq.(15) are our central results. Both of 
them have the form of a factorization law similar to 
Eq.(7) and Eq.(8). In Eq.(14), the factor is universal de- 
termined only by the channel's action on the maximally 
entangled state. But in Eq.(15), the factor includes -q rel- 
evant to the input state itself, which, however, is easy to 
compute by contrast to the evolution of the input state. 
So in order to know the dynamics of r of some pure in- 
put states, we only need to study the dynamics of r of 
the maximally entangled state and calculate the Schmidt 
coefficients of the input states, escaping from the cum- 
bersome task to compute the evolution equation of every 
different input state. Another fortunate thing is that 
unlike Eq.(8), the RHS of Eq.(14) and Eq.(15) can be 
calculated analytically. 

Here we would like to point out for a channel £ , if 
t((1 C§> £ )\4> + )(<fr + \) = 0, then for arbitrary input states, 
we simply find r((l ® £ )\ip)(ip\) = 0. However, we know 
(1 ® £)\ip)(ip\) may still be entangled since t is a lower 



bound of concurrence. In contrast to concurrence, if 
C[(l (g> £)\4> + )(<j) + \] — for a maximally entangled state, 
we know for arbitrary input states, C[(l (g> £)\ip) = 0, 
the output states are always separable. We know that £ 
is the entanglement breaking channel. 

We can generalize Eq.(14) to the case where the input 
state is a mixed state po- Suppose po has a decomposition 
that po — By the convexity of r, we have 

r((l <g> £)p ) = T{J2 iPi (l ® £)m(i/)i\) < EiJ>iT((l ® 

£M)m < fr((l®£)|4>+X4>+|)£ iPi C 2 (|^>). Con- 
sidering all decompositions of p 0l it's easy to see 

r((l ® £) Po ) < ^r((l ® 5)|0 + )(0+|)r'(p o ), (16) 

where t'(po) is the tangle of po and it has an easily com- 
putable formula for a bipartite mixed state po having no 
more than two nonzero eigenvalues [3] and some states 
with high symmetry like isotropic states 

In fact t' itself has a similar dynamical property to 
Eq.(14) and Eq.(16). In the following proof we neglect 
the normalization coefficients p, p' and p" for simplic- 
ity. First we suppose the input state is pure. If both 
p' and ps are pure states then according to the defini- 
tion of t' and Eq.(8) we have r'(p') < ^t' ( P£ )t' (\tp)). 
When p£ is mixed and has an optimal decomposition 

P£ = St-MViXVil such tnat T '(P£) = Si V^'G^i)); we 
have 

t'(p') <E a « t '((^® i)I^X^I) 

<jE A.r'd^yd^)) = t T ' {p£ y m . (17) 

i 

For the case of mixed input state, by the convexity of t', 
Eq.(17) also holds replacing \tp) with p a . 

V. EXAMPLES AND DISCUSSION 

Suppose £ is a depolarizing channel, such that £(p) = 
(1 — e)p + Egl with e e [0,1]- Using the definition 
in Eq.(4) to decompose (1 ® £ I mt° some 4x4 

matrices, we find that r((l ® £)|V')(V'|) = Si<j=o ^fi' 
where a;^ = max{0, 2d ~( 2 ^+ 2 ) e ^/uJiuTj } . Next we sup- 
pose £ is a phase damping channel, namely £{p) = 

(1 — e)p + e Y^,i=o Pu\i) {i\ for an input state p. Through 
calculation similar to that of the depolarizing channel, 

we obtain r((l <g> £)\ip) = Ei<i=Q3/y, wh ere = 
max{0, 2(1 — e) v /aJToJ7}. 

Now we study the case where the input state is mixed, 
for example an isotropic state pf = -felrj ( 1 — 1 4> + ) (0 + 1 ) + 
F\4>+)(4>+\, where F = (0+|p F |0+) e [0,1]. Due to 
its invariance under transformation T(p) = J dU(U ® 
U*)p{U ®U*)\ there exist elegant formulas for its tangle 
as well as concurrence [8] . Noting that if one qudit of pp 
goes through a depolarizing channel, pp is transformed 
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FIG. 1: The decay of r((l ® £)\ip)(ip\) (solid line) and its up- 
per(dashed line) and lower bound(dotted line), where (&,):£ is 
a depolarizing channel and (b) :£ is a phase damping channel. 
Here we let d = 3, u>o = ivi = ^ and = |. Note that in (a) 
a sudden death of r appears but in (b) it doesn't. 



FIG. 2: The decay of t((1 ® £)pF)(solid line) and its upper 
bound(dashed line), where £ is a depolarizing channel and 
d = 3. (a): F = §; (b): F = § . Note that r and its upper 
bound both vanish in finite time, although maybe not at the 
same time. 



into another isotropic state pf' with F' = F — Fd di x e 
and for isotropic states t(pf) is exactly equal with C 2 (/9p) 
[6], we have for F < ±, r((l (g) £ )p F ) = and for F > \, 

r[[l®£)p F ) = ^(max{0,F- \ - ^e}) 2 . 

We focus our attention on the dynamics of r (see Fig.l 
and Fig. 2). For depolarizing channel, we find that when 
e > g^j, r((l ® vanishes. A similar phe- 

nomenon appears for r((l(g)£ )/9f) when e > p%Zf - This 
is a sudden death of Tjjsimilar to the sudden death of en- 
tanglement Q [l(| [HI [13] ■ We hope r will not vanish in 
a finite time because then it can provide a non-trivial 
lower bound to squared concurrence and the state being 
evolving is still distillable [6]. 

Just like the sudden death of entanglement cannot ap- 
pear for any channel [9] , the sudden death of r doesn't 
exist for some channels. For example, for phase damping 
channel, we can see only when e = 1, r((l £g> £ )\tp)(ip\) 
vanishes, which means r doesn't die suddenly but asymp- 
totically. But, if the input state is mixed, for example a 
3® 3 Werner state, the sudden death of r can also appear 
even for phase damping channel. 

Summary. — The dynamics of a system is a fundamen- 
tal feature to describe its time evolution property. In 
this paper, we have shown the dynamical properties of 
the lower and upper bounds of squared concurrence re- 
spectively. Unlike the concurrence itself, the lower bound 
of the squared concurrence in this paper is computable. 
Thus our results are more reachable in various situations. 
We use depolarizing and phase damping channels as ex- 



amples and find r will vanish in finite time. Whether 
the entanglement sudden death appear depends both on 
the channel and the input state. Our result provides an 
easy way to estimate the dynamics of the entanglement 
in realistic physical systems. 
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